Application of a Self-Similar Pressure Profile to Sunyaev-Zel'dovich Effect 

Data from Galaxy Clusters 



Tony Mroczkowski,''^'^ Max Bonamente,'*-^ John E. Carlstrom,^-^'^-^ Thomas L. Culverhouse,^'' 

Christopher Greer,^'^ David Hawkins/° Ryan Hennessy,^-^ Marshall Joy,^ James W. Lamb/° 
Erik M. Leitch,^ '^ Michael Loh,^-'' Ben Maughan/^-'^ Daniel P. Marrone,^-^ " Amber Miller/ i'*''^ 
Stephen Muchovej,^ Daisuke Nagai/^'" Clem Pryke,*"-^ *^ Matthew Sharp,*"-^ and David Woody^° 

ABSTRACT 

We investigate the utility of a new, self-similar pressure profile for fitting 
Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters. Current SZ imaging 
instruments — such as the Sunyaev-Zel'dovich Array (SZA) — are capable of probing 
clusters over a large range in physical scale. A model is therefore required that can ac- 
curately describe a cluster's pressure profile over a broad range of radii, from the core 
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of the cluster out to a significant fraction of the virial radius. In the analysis presented 
here, we fit a radial pressure profile derived from simulations and detailed X-ray anal- 
ysis of relaxed clusters to SZA observations of three clusters with exceptionally high 
quality X-ray data: A1835, A1914, and CL 11226.9-^3332. From the joint analysis 
of the SZ and X-ray data, we derive physical properties such as gas mass, total mass, 
gas fraction and the intrinsic, integrated Compton j-parameter. We find that param- 
eters derived from the joint fit to the SZ and X-ray data agree well with a detailed, 
independent X-ray-only analysis of the same clusters. In particular, we find that, when 
combined with X-ray imaging data, this new pressure profile yields an independent 
electron radial temperature profile that is in good agreement with spectroscopic X-ray 
measurements. 

Subject headings: cosmology: observations — clusters: individual (Abell 1835, Abell 
1914, CL J1226.9+3332) — Sunyaev-Zel'dovich Effect 



1. Introduction 



The expansion history of the universe and the growth of large-scale structure are two of the 
most compelling topics in cosmology. As galaxy clusters are the largest collapsed objects in the 
universe, taking a Hubble time to form, their abundance and evolution are critically sensitive to 
the details of that expansion history. Cluster surveys can therefore provide fundamental clues to 



the nature and abundance of dark matter and dark energy (see, e.g. JWhite et al.ll 19931 : iFrenk et al. 
19991 : iHaiman et allboOll : Iweller et al.ll2002l) . 



While clusters have been extensively studied using X-ray observations of the hot gas that 
comprises the intraclus ter medium (ICM), radio me asurements of the same gas via the Sunyaev- 
Zel'dovich (SZ) effect (ISunvaev & Zerdovichlll972h provide an independent and complementary 
probe of the ICM (e.g., Carlstrom et al. 20021) . The SZ effect arises from inverse Compton scat- 



tering of cosmic microwave background (CMB) photons by the electrons in the ICM, imparting 
a detectable spectral signature to the CMB that is independent of the redshift of the cluster. As 
the SZ effect is a measure of the integrated line-of-sight electron density, weighted by tempera- 
ture, i.e., the integrated pressure (see § 12.11) . it probes different properties than the X-ray emission, 
which is proportional to the square of the electron density. Cluster surveys exploiting the redshift 
indep endence of the SZ effe ct are now being conducted by a variety of in struments, including the 
SZA (IMuchovei et all 120071). the South Pole Telescope dRuhl et all 120041) . the Atacama Cosmol- 
ogy Telescop e ( Kosowskyf 2003 ) and the Atacama Pathfinder Experiment (APEX) SZ instrument 
doobbs etal.ll2006 ). 
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To maximize the utility of clusters as cosmological probes we must understand how to ac- 
curately relate their observable properties to their total masses. The integrated SZ signal from a 
cluster is proportional to the total thermal energy of the cluster, and is therefore a measure of the 
underlying gravitational potential, and ultimately the dark matter content, within a cluster. The SZ 
flux of a cluster t hereby should provide a robust, low scatter proxy for the total cluster ma ss, Mjot 
(see, for example, Ida Silva et al.ll2004l : iMotl et allllOOsI : iNagailbood iReid & Spergellbooel) . 



To date, cosmological studies com bining SZA and X-ray data have relied a lmost exclusively 
on the isothermal /3-model (first used in Cavaliere & Fusco-Femianol 1976[ 1978) to fit the SZ sig 



nal in the region interior to r j^oo (see Grego et al. 200ol : 



Reese et al. 



2002 



LaRoque et al.ll2006l : 



Bonamente et al.ll2006L I2008L for applications of the isothermal ;S-model to SZ-i-X-ray data). Here 
^2500 is the radius within which the mean cluster density is a factor of 2500 over the critical density 
of the universe at the cluster redshift. While the isothermal j6 -model recovers global properties of 
clusters quite accurately in this regime ([LaRoque et al.ll2006l) . deep X-ray observations of nearby 
clusters show that isothermality is a poor description of th e cluster outskirts (r ~ rsoo) (see e.g. 
Piff'aretti et al.ll2005l: IVikhlinin et al.ll2005l : IPratt et al.ll2007L and references therein). An improved 
model for the cluster gas, accurate to large radii, is therefore critical for the analysis and cosmo- 
logical interpretation of SZ data obtained with new instruments that are capable of probing the 
outer regions (r ~ r^oo) of clusters. Such a model must be simple enough that it can be con- 
strained by SZ data with limited angular resolution and sensitivity typical of data sets acquired by 
SZ survey instruments optimized for detection, rather than imaging. While the jS-model has the 
virtue of simplicity, previous attempts to relax the assumption of isothermality typically required 
high- significance, spatially-resolved X-ray spectroscopy; such data are seldom obtained in short 
X-ra y exposures of high-r edshift clusters. This is particularly true for the cluster outskirts (see, 
e.g., LaRoque et al.ll2006|) . Attempts to move beyond the /3-model have typically improved the 
modeling of the gas density only within the core (r < 0.15 rsoo). 

In this work, we investigate a new model for the cluster gas pressure by using it to fit SZ data 
from the SZA and X-ray data from Chandra to thre e well-studied clust ers: Abell 1835, Abell 1914, 
and CL J1226.9-I-3332. The model was derived by iNagai et al.l (|2007|) from simulations and from 
detailed analysis of deep Chandra measurements of nearby relaxed clusters. The simplicity of 
this model — and the fact that SZ data are inherently sensitive to the integrated electron pressure — 
allow it to be used either in conjunction with X-ray imaging data, or fit to SZ data alone. The 
outline of the paper is as follows: in § [2l we present the details of the model, and couple it with 
an ICM density model that allows the inclusion of X-ray imaging data. In § [31 we apply this 
method to three clusters, combining new data from the SZA with Chandra X-ray imaging data. 
We demonstrate the utility of this model by applying it to SZ-i-X-ray data without relying on X-ray 
spectroscopic information. Results from the joint SZ-i-X-ray analysis are then compared to results 
from an X-ray-only analysis, including spectroscopic data, in §|4l We offer our conclusions in §[51 
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2. Cluster Gas Models 



2.1. Sunyaev Zel'dovich Effect and X-ray Emission 

The thermal SZ effect is a small (< 10"^) distortion in CMB intensity cause d by inverse Comp- 
ton scatterin g of CMB photons by energetic electrons in the hot intracluster gas (|Sunyaev & Zel'dovich 
197(aI1972|) . This spectral distortion can be expressed, for dimensionless frequency x = hvlk^Tcun, 
where h is Planck's constant, v is frequency, ks is Boltzmann's constant, and Tcmb is the primary 
CMB temperature, as the change A/^z relative to the primary CMB intensity normalization /q. 



A/ 
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Here ctj is the Thomson scattering cross-section of the electron, £ is the line of sight, and m^c^ 
is an electron's rest energy. The factor g(x, Tg) encapsulates the frequency dependence of the SZ 
effect intensity. For non-relativistic electrons. 
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We use the lltoh et al.l (|l998h relativistic corrections to Eq. [31 which are appropriate for fitting 
thermal SZ observations, and discuss their impact on the fit ICM profiles in § 14. 3[ Note that we 
have used the ideal gas law (Pe = neksTe) to obtain Eq. |2]from Eq. [B we use this to relate SZ 
intensity directly to the ICM electron pressure. 

The X-ray emission from massive clusters arises predominantly as thermal bremsstrahlung 
from electron s. The X-ray surface brightness produced by a cluster at redshift z is given by (e.g., 
Sarazin|[l988l) 

= TTT—Ta r '^l^eeiTe, Z) d£, (4) 

4;r(l -I- zy J 

where the integral is evaluated along the line of sight and the X-ray cooling function Ag^ is propor- 

1 /2 

tional to T,'\ The X -ray emission is therefore proportional to the square of the gas density, with 
a relatively weak dependence on T^. Separate spectroscopic observations of X-ray line emission 
can be used to measure the gas temperature. Thro ughout this work, we use the Raymond-Smith 
plasma emissivity code ([Raymond & Smithlll977b with constant metallicity Z = 0.3Zq, yielding 
He = 1-17 and /i = 0.61 for the mean m olecular weight of the elec trons and gas, respectively, using 
the elemental abundances provided by Anders & Grevesse ( 19891) . The SZ and X-ray observables 
A/5Z and S X are computed by evaluating the integrals in Eqs. [2] &|4] numerically. 
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2.2. Three-dimensional Models of the ICM Profiles 



In t his work, we adopt an analytic parameterization of the cluster radial pressure profile pro- 
posed bv lNagaietalJ (l2007h (hereafter N07), 



Pe{r) = 



ir/rpy[l+ir/rpy] 



{h-c)la ' 



(5) 



where Pe,i is a scalar normalization of the pressure profile, Vp is a scale radius (typically Vp « 
rjoo / 1 .3), and the parameters {a, b, c) respectively describe the slopes at intermediate (r « rp), outer 
(r > rp), and inner (r rp) radii. Note that Eq. [5]is a generalization of the analytic fitting for- 
mula obtain ed in numerical simu lations as a parameterization of the distribution of mass in a dark 
matter halo (|Navarro et al.lll997L NFW profile). This choice is reasonable because the gas pressure 
distribution is primarily determined by the dark matter potential. The use of a parameterized pres- 
sure profile is further motivated by the fact that self-similarity is best preserved for pressure, as 
demonstrated by the low cluster-to-cluster scatter seen when using these pa rameters to fit numeri- 
cal sim ulations. The NFW profile - in its pure form - has been applied by lAtrio-Barandela et al 



(|2008|) to fit the electron pressure profiles of SZ observations of clusters, who demonstrated an 
improvement over the application of the isothermal yS-model to SZ cluster studies. In this work, 
we adopt the fixed slopes (a, b, c) = (0.9, 5.0, 0.4), which NOV found to closely match the observed 
profiles of the Chandra X-ray clusters and the results of numerical simulations in the outskirts of 
a broad range of relaxed clusters I] 

The density model u s ed to fit the X-ray image data is a simplified version of the model em- 
ployed by Ivikhlinin et all (120061) (hereafter V06) is 



i.rlr,r 



1 



-I- 



^02 



[1 + {rlr,)^f~°'l^ [1 + {rlr,)yfly [1 + {rlr^^ff''' 



(6) 



We refer to Eq. |6l used in the independent X-ray analysis to which we compare our SZ-l-X-ray 
results, as the "V06 density model." Since the cluster core contributes negligibly to the SZ signal 
observed by the SZA, we exclude the inner 100 kpc of the cluster from the X-ray images used in the 
joint SZ-i-X-ray analysis. Recognizing a as the component introduced by IPratt & ArnaudI (|2002|) 
to fit the inner slope of a cuspy cluster density profile, and that the second /3-model component {^2) 
is present explicitly to fit the cluster core, we simplify the V06 density model to 



n,{r) = n,o [l + (r/r,)^]'""' [1 + (r/r,)^ 



-0.5e/y 



(7) 



'The original values published in N07 were {a, b, c) - (1 .3, 4.3, 0.7). These have recently been updated, however, 
and will be published in an erratum to N07. 
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where is the core radius, and is the radius at which the density profile steepens with respect to 
the traditional yS-model. Following V06, we fix y = 3, as this provides an adequate fit to all clusters 
in the V06 sample. We refer to this model as the "Simplified Vikhlinin Model" (SVM). Note that 
in the limit s ^ 0, Eq.|7]reduces to the standard jS-model. The SVM is thus a modification to the 
jS-model that has the additional freedom to extend from the core (r > r^) to the outer regions of 
cluster gas, spanning intermediate (r > r,) to large radii (r ~ r^oo). 

With the electron pressure and density in hand, we may also derive the electron temperature of 
the ICM using the ideal gas law, Te(r) = Pe{r)lkBne{r), where Pg and are given by Eqs.[5]and|71 
respectively. Note that Te{r) derived in this way is used in the analysis of X-ray surface brightness 
(Eq.lH). Hereafter, we refer to this jointly-fit cluster gas model as the N07-I-SVM profile. 



For comparison with previous work (e.g. lBonamente et al.ll2008h . we also employ the isother- 



mal jS-model for joint analysis of the SZ and X-ray data. In this model the density is given by 
Eq. |7] with £ = and Te{r) is a constant equal to the spectroscopically-measured temperature, 
Tx- The shape parameters of the isothermal ;0-model, Vc and jS, are jointly fit to the SZ and X-ray 
data, while the X-ray surface brightness (Eq. U) and SZ intensity profile (Eq. O normalizations are 
independently determined from the X-ray and SZ data, respectively. 



2.3. Parameter Estimation Using the Markov Chain Monte Carlo Method 

Our models have five free parameters to describe the radial distribution of the gas density (see 
Eq. |7]) and two parameters for the electron pressure (see Eq. [5]). Additional parameters such as 
the cluster centroid. X-ray background level, as well as the positions, fluxes and spectral indices 
of compact radio sources are also included where necessary. The Markov chain Monte Carlo 
( MCMC) method is used to extract the model parameters from the SZ and X-ray data, as described 



by iBonamente et al.l (120041) . In this section we provide a brief overview of this method, focusing 



on the changes to accommodate the NOV pressure model. 

The first step in fitting the SZ data is to compute the model image over a regular grid, sampled 
at less than half the smallest scale the SZA can probe. This image is multiplied by the primary 
beam of the SZA, transformed via EFT to Fourier space (where the data are naturally sampled by 
an interferometer; see § 14.11) . and interpolated to the Fourier-space coordinates of the SZ data. The 
likelihood function for the SZ data is then computed directly in the Fourier plane, where the noise 
properties of the interferometric data are well-characterized. 



The first step of the MCMC method is the calculation of the joint likelihood X of the X-ray 
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and SZ data with the model. The SZ likelihood is given by 



In(Xsz) = J] \~ (ARl + Alf) 



(8) 



where Ai?, and A/, are the differences between model and data for the real and imaginary compo- 
nents at each point / in the Fourier plane, and Wj is a measure of the Gaussian noise (1 /cr^). 

Since the X-ray counts, treated in image space, are distributed according to Poisson statistics, 
the likelihood of the model fit is given by 



H£nnage) = J] - M,- - ln(D, !)] 



(9) 



where M, is the model prediction (including cluster and background components), and Dj is the 
number of counts detected in pixel /. The inner 100 kpc of the X-ray images — as well as any 
detected X-ray point sources — are excluded from the fits by excluding these regions from the 
calculation of ln(X(mage) in Eq. HI 

The joint likelihood of the spatial and spectral models is given by X = -Csz -Cxmy For 
the N07-I-SVM fits, £xray is simply £image- Following Bonamente et all (120041 120061) . the X-ray 
likelihood for the yS-model fits is JLimage-Lxspec, as these must incorporate the likelihood Lxspec of 
the spectroscopic determination of Tx- The likelihood is used to generate the Markov parameter 
chains, and convergence of the chain to a stationary distribution i s established using the Raftery- 
Lewis and Geweke tests (IRafterv & Lewisll 19921 : bilks et al.lll996l) . 



2.4. Calculation of Mgas, Mtot, and Tint 

With knowledge of the three-dimensional gas profiles, we compute global properties of galaxy 
clusters as follows. The gas mass Mgas(r) enclosed within radius r is obtained by integrating the 
gas density Pgas = figinprieir) over a spherical volume: 



Mgas(r) = 4n[ p^,,{r'y^dr' . (10) 
Jo 

The total mass, Mtot, can be obtained by solving the hydrostatic equilibrium equation as: 

r2 dP{r) 

MUr) = -- r--r^, (11) 

where P = (jj.e/fj)Pe is the total gas pressure. We then compute the gas mass fractions as /gas = 

Mgas /Mtot. 
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The line of sight Compton j-parameter, which characterizes the strength of Compton scatter 
ing by electrons, is defined 

kBCTT 



y = 



■/"■ 



(12) 



mgC^ J 

We compute the volume-integrated Compton j-parameter, Y, from the pressure profile fit to the SZ 
observations for both cylindrical and spherical volumes of integration. The cylindrically-integrated 
quantity, Fcyi, is calculated within an angle 6 on the sky, corresponding to physical radius R = 6dA 
at the redshift of the cluster, 

y{e) ffdff = In I y{r') r'dr' = ^ I d€ I P,{r') r'dr'. (13) 

Jo '^eC J_oo Jo 



The last form makes explicit the infinite limits of integration in the line of sight direction, originat- 
ing with the definition of y (Eqs.[T]and[2l). The spherically-integrated quantity, Ygph, is obtained by 

rwr\fi\c^ \\r^i^r\^r\ q forliiic v ft*r\m i)r\c^ r»liictf:»t* nc^nif^-r 



integrating the pressure profile within a radius r from the cluster center. 



Jo 



Pe{r') r"dr. 



(14) 



The parameter ysph('') is thus proportional to the thermal energy content of the ICM. 



To compute the global cluster properties descri bed above, one needs to d efine a radius out 



to whi ch all quantities will be calculated. Following iLaRoque et al.l (|2006|) and iBonamente et al 



(120061) . we compute global properties of clusters enclosed within the overdensity radius r^, within 
which the average density (p) of the cluster is a specified fraction A of the critical density, via 



-npc(z) Ari = Mtot(rA), 

where pdz) is the critical density at cluster redshift z, and A = 
evaluate cluster properties at density contrasts of A = 2500 and A 
r25on has been used in previou s OVRO and BIMA interferometric SZ studi es (e.g. iLaRoQue et al 



(15) 

{p)/Pc(z)- In this work, we 
= 500. The ov erdensity radius 



Allen et al 



20061: IBonamente et al.ll2008h as well as in many X-ray cluster studies (e.g. IVikhlinin et al.ll2006l : 



2007b . while rjoo is a radius reachable with SZA and deep Chandra X-ray data. 



3. Data 



3.1. Cluster Sample 



For this work, we selected three massive clusters that are well studied at X-ray wavelengths, 
and span a range of redshifts (z = 0.17-0.89) and morphologies, on which to test the joint analysis 
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of the Chandra X-ray and SZA data. We assume a ACDM cosmology with Qm 
and h = 0.7 throughout our analysis. 



0.3, Qa = 0.7, 



Located at z = 0.25, Abell 1835 (A1835) is an intermediate-redshift, rela xed cluster, as 



evide nced by its circular morphology in the X-ray images and its cooling core (|Peterson et al. 



200 lb . To demonstrate the applicability of our technique for high redshift clusters, we analyzed 
CL J 1226.9-H3332 (CL1226), an apparently relaxed cluster at z = 0.89 JMaughan et all 120041 
20071) . To assess how this method performs on somewhat disturbed clusters, we also analyzed 
Abell 1914 (A1914), an intermediate redshift {z = 0.17) cluste r with a hot subclump near the 
core. When the subclump is not excluded from the X-ray analysis. iMaughan et al.l(|2008l . hereafter 



M08) find a large X-ray centroid shift in the density profile, which they use as an indicator of an 
unrelaxed dynamical state. 

In the following sections, we discuss the instruments, data reduction, and analysis of the SZ 
and X-ray data. Details of these observations, including the X-ray fitting regions, the unfiagged, 
on-source integration times, and the pointing centers used for the SZ observations, are presented 
in Tables [U and [21 We also review an independent, detailed X-ray-only analysis, with which we 
compare the results of our joint SZ-i-X-ray analyses. 



3.2. Sunyaev-Zel'dovich Array Observations 

The Sunyaev-Zel'dovich Array is an interferometric array comprising eight 3.5-meter tele- 
scopes, and is located at the Owens Valley Radio Observatory. For the observations presented 
here, the instrument was configured to operate in an 8-GHz-wide band covering 27-35 GHz using 
the 26-36 GHz receivers (hereafter referred to as the "3 0-GHz" band) or cover ing 90-98 GHz us- 



ing the 80-1 15 GHz receivers (the "90-GHz" band). See lMuchovej et all (|2007l) and Marrone et al. 
(in preparation) respectively for more details about commissioning observations performed with 
the 30-GHz and the 90-GHz SZA instruments. Details of the observations presented here, includ- 
ing the sensitivity and effective resolution (the synthesized beam) of the long and short baselines, 
are given in Table [T] 

An interferometer directly measures the amplitude and phase of Fourier modes of the sky 
intensity, with sensitivity to a range of angular scales on the sky given by ~ A/B, where B is the 
projected separation of any pair of telescopes, i.e., a baseline. The field of view of the SZA is given 
by the primary beam of a single telescope, approximately 10.7' at the center of the 30-GHz band. 
At 30 GHz, optimal detection of the arcminute-scale bulk SZ signal from clusters requires the short 
baselines of a close-packed array configuration; six of the SZA telescopes were arranged in this 
configuration for the observations presented here, yielding 15 baselines with sensitivity to ~ 1-5' 



scales. The two outer antennas, identical to the inner six, provide an additional 13 long baselines in 
this configuration, with sensitivity to small-scale structure, allowing simultaneous measurement of 
compact radio sources unresolved by the long baselines (angular size < 20") which could otherwise 
contaminate the SZ signal. Observations at 90 GHz with the SZA probe scales at three times the 
resolution of the 30-GHz observations for the same array configuration. The short baselines of the 
90-GHz observations thereby bridge the gap between long and short baseline coverage at 30-GHz. 

SZA data are processed in a complete pipeline for the reduction and calibration of interfero- 
metric data, developed within the SZA collaborati on. Absolute flux calibrations are derived from 
observations of Mars, scaled to the predictions of iRudyl (Il987|) . Data from each observation are 
bandpass-calibrated using a bright, unresolved, flat-spectrum radio source, observed at the start or 
end of an observation. The data are regularly phase-calibrated using radio sources near the targeted 
field; these calibrators are also used to track small variations in the antenna gains. Data are flagged 
for corruption due to bad weather, sources of radio interference and other instrumental effects that 
co uld impact data quality . A more detailed account of the SZA data reduction pipeline is presented 



in 



Muchovej et all (l2007h . 



In the SZA cluster observations presented here, the A1835 field contains three detectable 
compact sources at 31 GHz: a 2.8 + 0.3 mJy central source, a 1.1 + 0.4 mJy source about one 
arcminute from the cluster center, and a 0.8 + 0.4 mJy source 5.5 arcminutes from center. The po- 
sitions of these sources are in good agreement with those from the NVSS (which only contains the 
central source) and the FIRST surveys. The SZA observation of CL1226 contains one detectable 
compact source, identified in both FIRST and NVSS, wi th flux at 31 GHz of 3 .9 + 0.2 mJy. This 
source is 4.3 arcminutes from the cluster center (see also iMuchovej et al.ll2007h and therefore lies 
outside the field of view of the SZA 90 GHz observations. Two compact sources, with positions 
constrained by NVSS and FIRST, were detected at 31 GHz in the A1914 field. The fluxes of these 
sources are 1.8 ± 0.4 mJy and 1.2 + 0.3 mJy. The stronger was detected in both the NVSS and the 
FIRST surveys, while the weaker was only detected in the FIRST survey. Tabled] summarizes the 
properties of the compact radio sources detected in the SZA observations. 

When fitting compact sources detected in the SZ observations, we calculated the spectral 
indices fror n the measured flux a t 3 1 GHz a nd 1.4 GHz, where the latter was constrained by either 
the NVSS dCondon et all 1 19981) or FIRST dwhite et al.lll997l) survey, respectively. The source 
fluxes an d approximate c oordinates are first identified using the interferometric imaging package 
Difmap (IShepherdlll997h . We first refine the source positions by fitting a model in the MCMC 
routine, and then fix these positions to their best-fit values when fitting the cluster SZ model. We 
leave the source flux a free parameter, so that the cluster SZ flux and any compact sources are fit 
simultaneously. 
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3.3. X-ray Observations 

All X-ray imaging data used in this analysis were obtained with the Chandra ACIS-I detector, 
which provides spatially resolved X-ray spectroscopy and imaging with an angular resolution of 
~ 0.5" and energy resolution of ~ 100-200 eV. Table [2l summarizes the X-ray observations of 
individual clusters. 

For the X-ray data used in the joint SZ-l-X-ray analysis, images were limited to 0.7-7 keV in 
order to exclude the data most strongly affected by background and by calibration uncertainties. 
The X-ray images — which primarily constrain the ICM density profiles — were binned in 1 .97" 
pixels; this binning sets the limiting angular resolution of our processed X-ray data, as the Chandra 
point response function in the center of the X-ray image is smaller than our adopted pixel size. The 
X-ray background was measured for each cluster exposure, using source-free, peripheral regions 
of the adjac ent detector ACIS - I chips. Ad ditional details of the Chandra X-ray data analysis are 



presented in lBonamente et al.l (12004 120061) . 



lip: 

3, 



In § 14.31 & 14.41 we compare the results of our joint SZ-l-X-ray analysis to the results of inde- 
pendent X-ray analyses. For A 1914 and CL1226 we use the data and analysis described in detail in 
M08 and lMaughan et all (120071) (hereafter M07), respectively. For A1835, the ACIS-I observations 
used became public after the M08 sample was published, and we therefore present its analysis here 
for the first time. The observation of A1835 was calibrated and analyzed using the same methods 
and routines described in M08, which we now briefly review. 

In the X-ray-only analyses, blank-sky fields are used to estimate the background for both the 
imaging and spectral analysis. The imaging analysis (primarily used to obtain the gas emissivity 
profile) is performed in the 0.7-2 keV energy band to maximize signal to noise. Similar to the joint 
SZ-l-X-ray analyses, these images were also binned into 1.97" pixels. 

For the spectral analysis, spectra extracted from a region of interest were fit in th e 0.6-9 keV 



band with an absorbed, redshifted Astrophysical Plasma Emission Code (APEC) dSmith et al 



200 1|) model. This absorption was fixed at the Galactic value. This spectral analysis was used 
to derive both the global temperature Tx, determined within the annulus r G [0.15, 1.0] r^oo, and 
the radial temperature profile fits of the V06 temperature profile, given by 



(r/r^ooi)"^'"" + 1 



{r/r,r 



(1 + (r/rr)'^) 



by/h 



(16) 



We refer the reader to V06 for details, but note that this temperature profile is the combination of a 
cool core component (first set of square brackets, where the core temperature declines to r^in) and 
a decline at large radii (second set of square brackets, where temperature falls at r > r^). 



An important consideration when using a blank-sky background method is that the count rate 
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at soft energies can be significantly different in the blank- sky fields than in the target field, due 
to differences between the level of the soft Galactic foreground emission in the target field and 
that in the blank-sky field. This was accounted for in the imaging analysis by normalizing the 
background image to the count rate in the target image in regions far from the cluster center. In 
the spectral analysis, this was modeled by an additional thermal component that was fit to a soft 
residual spectrum (the di fference between spec tra extracted in source free regions of the target and 



background datasets; see lVikhlinin et al.l (|2005l) ). The exception to this was the XMM-Newton data 
used in addition to the Chandra data for CL1226. As discussed in MOV, a local background was 
found to be more reliable for the spectral analysis in this case, thus requiring no correction for the 
soft Galactic foreground. 

The M07/M08 X-ray analysis methods exploit the full V06 density and temperature models 
(Eqs.[6]&[T6l respectively) to fit the emissivity and temperature profiles derived for each cluster, 
and the results of these fits are used to derive the total hydrostatic mass profiles of each system. 
Uncertainties for the independent. X-ray-only analysis method are derived using a Monte Carlo 
randomization process. These fits involved typically ~1000 realizations of the temperature and 
surface brightness profiles, fit to data randomized according to the measured noise. We refer the 
reader to M07 and M08, where this fitting procedure is described in detail. 



4. Results 
4.1. SZ Cluster Visibility Fits 



Interferometric SZ data are in the form of visibilities Vy{u, v) (see, for example. [Thompson et al 



200 ih . which for single, targeted cluster observations with the SZA can be expressed in the small 



angle approximation as 

V,{u, v) = JJ Mx, y) lyix, y) e-'2-(-+^3') dx dy. (17) 

Here u and v (in number of wavelengths) are the Fourier conjugates of the direction cosines x and y 
(relative to the observing direction), Ay{x, y) is the angular power sensitivity pattern of each antenna 
at frequency v, and Iv{x,y) is the intensity pattern of the sky (also at v). Eq.[l7]is recognizable as a 
2-D Fourier transform, so the visibilities give the flux for the Fourier mode for the corresponding 
u, v-coordinate. 

By combining Eqs. [Hand [171 we can remove the frequency dependence from the measured 
cluster visibilities, just as we have related SZ intensity to the frequency-independent Compton 
j-parameter. We define the frequency-independent cluster visibilities Y(u, v) as 

Vy(u,v) = gix)IoY(u,v), (18) 
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where Iq (in units of flux per solid angle) is 

{hcY 

Additionally, we rescale Y{u, v) by the square of the angular diameter distance, J^, in order to 
remove the redshift dependence from the cluster SZ signal. Note that, while we use the non- 
relativistic g(x) (Eq. O to compute Y{u,v) (Eq. [TS]) , we only use this for display purposes. The 
effects of assuming the classical SZ frequency dependence are discussed in § 14.31 

Figure [H shows the maximum-likelihood fits of the N07 profile and the isothermal yS-model to 
each cluster's visibility data, from which we have subtracted the detected radio sources (Table |3]). 
We also removed the frequency dependence of the SZ eff"ect by rescaling the cluster visibilities to 
Y(u, v) (Eq. [T8l). For the purposes of plotting, this rescaling is useful when binning the SZ signal 
across 8 GHz of bandwidth as well as when plotting the 30-GHz and 90-GHz SZA data taken on 
CL1226. 

As indicated in Fig. [H both the isothermal yS-model and NOV model (which was fit jointly 
with the SVM) fit the available data equally well. However, as Vw^ + (large scales on the 
sky), the isothermal jS-model extrapolates to a much larger value of Y(u,v). This corresponds to 
the much larger values of Fcyi that are computed at large radii using fit jS-models (see §14.41) . 

The middle panel of Fig. [T] shows the combined 30+90 GHz observations of CL1226, which 
has a smaller angular extent than A 1835 or A1914 due to its distance (compare, for example, the 
values of r^oQ for each cluster, listed in Table |5]). The dynamic range and M,v-space coverage pro- 
vided by the 30-GHz SZA observations (black points) on CL1226 were insufficient for constrain- 
ing the radial pressure profile of the cluster. The short baselines of the 90-GHz SZA observations 
(middle three points) help to provide more complete m,v- space coverage, as discussed in ^3.2[ 



4.2. X-ray Surface Brightness Fits 

The X-ray surface brightness (Eq.®, ignoring the data within a 100 kpc radius, was modeled 
separately with both the isothermal yS-model, using the spectroscopically-determined, global Tx 
(measured within r e [0.15, 1.0] rsoo), and the SVM, using the temperature derived from the NOV 
pressure profile fit to the SZ data. Figure [2] shows the maximum-likelihood fits to the surface 
brightness of each cluster for both the SVM and isothermal ;0-model. For plotting purposes, the X- 
ray data are radially-averaged around the cluster centroid, which is determined in the joint SZ-i-X- 
ray analysis by fitting the two-dimensional X-ray imaging data with the spherically-symmetric 
SVM and isothermal yS-model profiles. 
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4.3. ICM Profiles 

Figure [3] shows the three-dimensional ICM radial profiles derived from the joint analysis of 
SZA + Chandra X-ray surface brightness data for A1835, CL1226 and A1914 (from left to right). 
From top to bottom, we show the electron pressure, the gas density and the derived electron tem- 
perature profiles, each as a function of cluster radius. In all panels, we compare the ICM profiles 
derived from the N07-I-SVM model to the results of a traditional isothermal yS-model analysis, indi- 
cated by solid and dot-dashed lines, respectively. The hatched regions indicate the 68% confidence 
interval for each derived parameter. 

As shown in the top panels of Fig. [3l the pressure profiles derived from the NOV model and 
the isothermal yS-model show good agreement within r25oo, but deviate by ~3-5-cr in the cluster 
outskirts. This is a consequence of the fact that clusters exhibit a significant decline in tempera- 
ture beyond r25oo, as determined from spectroscopic X-ray observations (see the bottom panel of 
Fig. [3]). The pressure profile derived from the isothermal jS-model analysis is therefore biased sys- 
tematically high beyond r25oo- In contrast, the NOV model, which fits the pressure directly, is free 
to capture the true shape of the pressure profile well beyond the radius at which the assumption of 
isothermality becomes invalid. 

For all three clusters there is little evidence for the second component of the electron density 
allowed by the SVM; the density is fit equally well by either the SVM or a single-component /3- 
model, as illustrated by the center row of panels in Fig. [3] (see also Table [2l). For all three clusters, 
the fits of the SVM agree to within 1-2% of the full V06 density profile fits (not shown) outside 
the core; discrepancies at this level are easily attributed to differences in the fitting of the X-ray 
background, and to the differences between the APEC and Raymond-Smith emissivity models 
used respectively in the M0V/M08 X-ray analysis and the joint SZ-l-X-ray analyses. 

In the bottom panels, the electron temperature profiles inferred from the NOV-i-SVM pro- 
files are compared to temperature profiles derived from deep spectroscopic Chandra X-ray ob- 
servations (and XMM-Newton in the case of CL1226; see MOV). The comparison shows that the 
radial electron temperature profiles derived from the NOV-i-SVM profiles are in good agreement 
with independent X-ray measurements for the relaxed clusters A1835 and CL 1226, which exhibit 



radially decreasing te mperature profiles in the cluster outskirts (see also iMarkevitch et al.lll998 



Vikhlinin et al.l 120051) . The disturbed cluster A1914, however, shows less overall agreement be- 
tween the derived NOV-l-SVM radial temperature profile and the M08 fit of the V06 temperature 
profile. Since the NOV pressure profile fit to the SZA observation of A1914 agrees with that derived 
from M08 within their respective 68% confidence intervals, the temperature discrepancy is likely 
due to deviations from the spherical symmetry implicitly assumed in this analysis. Additionally, 
scales greater than about six arcminutes are beyond the radial extents probed by the SZA; it is 
unsurprising the agreement becomes poorer at radii larger than this. 
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We note that we fit the SZ da ta using relativistic corrections to the SZ frequency depen- 
dence provided by lltoh et alj (|1998b . These corrections are appropriate for the thermal SZ effect 
at the SZA observing frequencies of 30 and 90 GHz. Compared to fits assuming the classical SZ 
frequency dependence, a pressure profile fit using the relativistic corrections has both a higher nor- 
malization and larger upper error bars. This is noticeable when including higher frequency data, 
where the relativistic correction is larger (~ 5% at 90 GHz versus ~ 3% at 30 GHz for cluster tem- 
peratures ~ 8 keV). This increase in the pressure fit is due to the diminished magnitude of the SZ 
effect w hen using t hese c orrections (for frequencies below the null in the SZ spectrum, < 218 GHz; 
see e.g. lltoh et al.l (|l998|) ). The pressure profile therefore must adjust to fit the observed SZ flux. 



The larger upper error bars on the fit pressure profile arise from a more subtle effect. Since the 
temperature is derived from the simultaneously fit pressure and density profiles, and the SZ effect 
diminishes as electrons become more relativistic (i.e. hotter), the upper error bar of the pressure 
fit must increase to fit the same noise in the observation (compared to the non-relativistic case). 
The lower error bar is less affected, as lower electron temperatures require smaller relativistic 
corrections. The resulting asymmetric error bars can be seen in the derived temperature profile of 
CL1226, which relied on 90 GHz data, in Figured 



4.4. Measurements of Y, Mgas, Mtot, and /gas 

In Tables |l]and[5l we report global properties of individual clusters derived from the N07-I-SVM 
model fits to the SZ-l-X-ray data. We calculate all quantities at overdensity radii r25oo and rjoo, and 
compare them to results from both the isothermal /3-model analysis of the same data, as well as to 
the X-ray-only analysis. 

The NOV pressure profile has just two free parameters — Pej and — which exhibit a degener- 
acy. Figure |4] shows this degeneracy in fits of the N0 7 profile to the SZA observations of A1835. 



Similar to the -/3 degeneracy of the jS-model (see iGrego et al.ll200lL for example), these two 
quantities are not individually constrained by our SZ observations, but they are tightly correlated 
and the preferred region in the Pej - rp plane encloses approximately constant Ycyi- As a result, 
the 68% confidence region for Y^yi is more tightly constrained than the large variation in P^ i or 
Tp might individually indicate. Figure |4] also shows that the inclusion of X-ray data has only a 
marginal effect on the value of Fcyi derived from the SZ fit. This is as expected, due to the weak 
dependence of the X-ray surface brightness on temperature (see § 12.11 ) and the fact that the NOV 
profile is not linked to the SVM density profile's shape. This indicates that X-ray data are not 
necessary to constrain Y^yu although they do limit the range of accepted radial pressure profiles. 

At both r25oo and rjoo, the measurements of Ycyi derived from the joint N07-I-SVM and the 
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X-ray-only analysis are consistent at the 1-cr level, for all three clusters. The isothermal ;8-model 
analysis, however, overestimates Ycyi by ~ 20%-40% at r25oo, and by ~ 30%-l 15% at rjoo- This is 
due to the large contribution to Ycyi from the cluster outskirts, at every projected radius, where the 
jS-model significantly overestimates the pressure. 

In contrast to the systematic variations in Ycyh the determinations of Fsph(r25oo) are consistent 
among the three analyses. At r^Qo, however, the median Y^p^ values from the isothermal /3-model 
can be as much as ~ 60% higher than either the N07+SVM or M08 results, due to the fact that 
isothermality is a poor description of the cluster outskirts. 

The gas mass estimates computed using either the jointly-fit SVM or the isothermal yS-model 
agree with the gas masses derived from the Maughan X-ray fits (Tables |4] and [5]). This agreement 
is not surprising, given that the gas mass is determined in all cases from density fits to the X- 
ray data. It demonstrates, however, that the 100 kpc core makes a negligible contribution to the 
total gas mass even at r25oo, and that excluding the core from the joint analysis does not therefore 
introduce any significant bias in our estimate of Mgas. Incidentally, it also shows that the additional 
component allowed by the SVM and the full V06 density models is not indicated in these clusters. 

In Tables |4] and [5l we also present estimates of the total masses, computed using each model's 
estimate of the overdensity radius (r/^) for each Monte Carlo realization of the fit parameters. For 
two of the clusters, we find that the error bars are significantly larger for Mtot determined from the 
N07-I-SVM fits than for the isothermal /3-model or M08 fits. This is a consequence of the fact that 
the yS-model analysis, with fewer free parameters and the assumption of isothermality, eff"ectively 
places stronger but poorly-justified priors on the total mass. We find that the N07-I-SVM and M08 
total mass estimates agree at both rsoo and r25oo, leading to good overall agreement between gas 
fractions computed using the NOV-i-SVM profiles and those from the Maughan X-ray fits. 

Since the HSE estimate for Mtot is sensitive to the change in slope in the pressure profile 
(dP/dr), Mtot is not as well constrained by the N07-I-SVM profiles as is Y^yu which scales directly 
with integrated SZ flux, a parameter that is more directly linked to what the SZA observes (see 
§ 14.11) . Figure |5] shows a comparison of the NOV-i-SVM estimates for Mtot and Fcyi, revealing that 
Ycyi has a more tightly constrained and centrally peaked distribution than Mtot does. 

The isothermal /3-model, on the other hand, is over-constrained such that it cannot agree with 
the non-isothermal estimates of Mtot at both r25oo and r^oo. Its estimates of Mtot are moreover sen- 
sitive to the annulus within which Tx is determined (see § 13.31 ). This trend can be seen in Fig. [6l 
which shows Mtot(?") for each cluster. 
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5. Conclusions 



In this work, we demonstrated the application of a new pressure profile, proposed by lNagai et al 



(l2007h . in fitting SZ data taken by the Sunyaev-Zel'dovich Array. By combining the pressure pro- 



file constraints from the SZ data with constraints on density from Chandra X-ray imaging, using 



a simplified form of the density model proposed by iVikhlinin et alj (120061) . we were able to de- 
termine the properties of three clusters (A1835, A1914, and CL J 1226.9-1-3332) spanning a wide 
range in redshift and dynamical states. 

This technique was compared with two others: a joint analysis of the same SZ -I- X-ray data 
relying on the commonly used isothermal jS-model, which assumes the pressure and density profiles 
have the same form, and an X-ray only analysis that independently models temperature and density 
using both X-ray imaging and spectroscopic measurements (but excluding SZ data). We find that 
the global cluster properties at both r25oo and rjoo determined from ICM profiles fit in the joint 
pressure analysis are in excellent agreement with properties determined in the X-ray analysis. By 
contrast, the isothermal yS-model tends to overestimate with respect to the other models the gas 
pressure at r > rasoo, where isothermality is an increasingly poor assumption. The yS-model thus 
leads to an overestimate of the cylindrically-integrated Compton y-parameter at rsoo- Since the 
isothermal ;0-model does not provide a good description of the ICM profile in the cluster outskirts, 
we caution against its use in deriving global properties of clusters even at r ~ rjoo, which is a large 
fraction of the virial radius. 

We tested the ability to recover the ICM electron pressure profiles from SZ data by analyzing 
the SZ data together with the X-ray imaging data alone, ignoring the X-ray spectroscopic informa- 
tion. Assuming the ideal gas law, we derive electron temperature profiles by coupling the pressure 
fits to the SZ data with the density fits to the X-ray imaging data. We find that these derived 
temperature profiles show broad agreement with those determined spectroscopically from deep 
X-ray observations, even for the highest redshift cluster in our sample, at z = 0.89. This method 
therefore provides an independent technique for determining the radial electron temperature distri- 
bution in high-redshift clusters, for which deep spectroscopic X-ray data may be unavailable and 
are impractical to obtain. 
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Table 1 . SZA Cluster Observations 



Cluster Name 


2" 


Da 


SZA Pointing Center (J2000) 


Short Baselines (0.3-1.5 U) 


Long Baselines (3.0-7.5 k/i) 


^int,red 






(Gpc) 


a 


6 


beam(" X "f 


o-(mJy)'' 


beam(" X "f 


o-(mJy)'' 


(hi-s) 


A1914 


0.17 


0.60 


14''26'"00'.8 


+37°49'35".7 


117.5x129.9 


0.30 


23.5x17.4 


0.35 


11.5 


A1835 


0.25 


0.81 


14''01'"02'.0 


+02°52'41".7 


116.6x152.1 


0.25 


17.5x23.5 


0.33 


18.6 


CL1226 (30-GHz) 


0.89 


1.60 


12''26"'58'.0 


+33°32'45".0 


117.4x125.4 


0.20 


16.0x21.2 


0.20 


22.0 


CL1226 (90-GHzf 










42.3x39.1 


0.42 






29.2 



"Redshifts for A1914 and A1835 are from lStruble & RoodI 1 19991) . Redshift for CL1226 is from lEbeling et alj bOOlh . All are in agreement with 
XSPEC fits to iron emission lines presented in lLaRoqueeT^ 20061) . 

''Unflagged data after reduction. 

"^Synthesized beam FWHM and position angle measured from North through East 
''Achieved rms noise rn coiTesponding maps 

■^The short baselines of the 90-GHz observation probe ~ 1—4.5 k/i. The long baseUnes of the 90-GHz observation were not used here. 



Table 2. Details of X-ray Observations 



Cluster Name 


ObsID 


Joint SZ+X-ray Analysis 


Maughan X-ray Analysis 






Tint 


fitting region'' 


' int 


fitting region 






(ks) 


ff 


(ks) 




A1914 


542+3593 


26.0 


34.4-423.1 


23.3 


0.0-462.5 


A1835 


6880 


85.7 


25.6-344.4 


85.7 


0.0-519.6 


CL1226 


3180+5014 


64.4 


12.9-125.0 


50 


0.0-125.0 




0200340101 


N/A 


N/A 


75+68"= 


17.1-115"= 



^Good (unflagged, cleaned) times for X-ray observations after respective calibration pipelines. 

''X-ray image fitting region. Excluded regions are not used in the X-ray likelihood calculation (Eq.|9]l. 

'^XMM-Newton observation of CL1226, presented in M07, was used only in the spectroscopic analysis. The exposure 
times are, respectively, those of the MOS and PN camera. 



-23- 



Table 3. Unresolved Radio Sources in 30-GHz Observations 



Cluster Name Src # Differential Offset from Pointing Center Flux at 30.938 GHz 

mJy 



A R.A. (") 



A Dec. (") 



A1914 


1 


-242.1 


-235.4 


1.8 


A1914 


2 


-160.5 


-108.6 


1.2 


A1835 


1 


0.9 


1.5 


2.8 


A1835 


2 


-22.8 


-51.2 


1.1 


A1835 


3 


-275.1 


-178.6 


0.8 


CL1226^ 


1 


258.4 


-38.9 


3.9 



+0.4 
-0.4 
+0.3 
-0.3 
+0.3 
-0.3 
+0.4 
-0.3 
+0.4 
-0.4 
+0.2 
-0.2 



''No point sources were detected in the 90-GHz observation of CL1226. 



Table 4. Y, 



cyl: 



t, and /gas for Each Model Computed Within r25oo- 



Cluster 


^2500 


'"2500 


icyl 


Ysph 




/gas 


model fit 


(") 


(Mpc) 


(10-^Mpc-) 


(lO-^Mpc^) 


(lO'^Mo) (lO'^Mo) 





AbeU 1835 

N07-HSVM 
Maughan (this work) 
isothermal yS-model 

CL J1226+3332.9 

N07-HSVM 
Maughan et al. (2007) 
isothermal jS-model 

AbeU 1914 

N07+SVM 
Maughan et al. (2008) 
isothermal yS-model 



173+5.8 
^ '•'-5.7 

169+5.5 
159+3.0 



52.8 



-1.8 
-1.9 

57 3^'-^ 



228 



-12.9 
-11.7 

218!^:' 

+5.7 



204 



0.68!°:°^ 
0.66t- 



41+0.01 
^■^^-0.01 

45+0°' 
"■^-^-0.01 

42+004 

"-'■^^-0.03 



67+0 0'* 

-0.03 

63+002 

"■"■'-0.02 

+0.02 



0.59 



11.73![i^ 
13.85!":^? 



c 49+0.53 
•'■^^-0.50 
7 C7+O.33 
'■■"-0.34 

6 76+0-''i 

"■'"-0.73 



8 43+'-^^ 

''■^•'-1.24 

7.87t^^ 
11.28-- 



o 9C+O.8I 
"■^-'-0.78 

7.88!0;4^ 
7 Q4+0.43 



3 56+0-^'' 

■'■•'"-0.36 

5 04+0-3 1 

■^■"^-0.28 
T, Q9+0.60 
■-'■='^-0.48 



+ 1.03 

-0.82 
a+0.37 
-0.38 

6 24+0-34 

"■^^-0.32 



6.29! 
5.69^ 



c 74+0.26 
''■ '^-0.25 
c 77+0.25 

J. / /_Q 35 

4 96+0- '3 
^■^"-0.12 



9 Qi +0.14 
^■^^-0.14 
a 9C+O.14 
■^■■^-0.13 

9 QT+0.40 

^■^■-^~0.33 



4 71+0.34 
^- ' ^-0.30 

4 64+0 " 

^■" -0.16 
+0.15 



4.04 



5 64+0-58 

-"■"^-0.54 

5 30+0-" 
4 38+0-25 

-0.24 



2 67+0-29 

-0.27 
T. 41+0.30 
■^■^^-0.26 
9 oQ+0.90 
^■°^-0.63 



4 97+0.89 
^■^'-0.72 
4 31+0.43 

t.Ji _Q 33 

o c 9 +0.30 
■'■■'^-0.26 



102+0006 

"■^"^-0.005 

109+0009 

"■^"^-0.006 

n 1 1 T +0.004 

"■^ ^■'-0.004 



109+0007 

"■^"^-0.007 

095+0-004 

"■"^•'-0.004 

101+0-014 

"■^"^-0.014 



095+0-009 

"■"^-'-0.009 

107+0-005 

"■^"'-0.006 

115+0-005 

"■^^•'-0.005 
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Table 5. Ycyh Y^ph, Mg^s, M^ot, and /gas for Each Model Computed Within r^oQ. 



Cluster 
model fit 



^500 
(") 



'"500 

(Mpc) 



(lO-^Mpc^) 



i'sph 

(lO-^Mpc^) 



(lO'^Mo) 



Mtot 
(lO'^'Mo) 



gas 



Abell 1835 

N07+SVM 
Maughan (this work) 
isothermal yS-model 

CL J1226+3332.9 

N07+SVM 
Maughan et al. (2007) 
isothermal yS-model 

Abell 1914 

N07+SVM 
Maughan et al. (2008) 
isothermal ;8-model 



369 



+28 
27 



361!^ 



-13 
-10 

115+^ 
hll 

-9 



127 



126 



430!^^ 
461!|3 



1 44+0.11 
^■^^-0.10 
1 49+0.07 
^ -^^-0.05 
1 41+0.03 
^■^^-0.03 



98+0.10 

^■^°-0.07 

+0.02 
-0.02 
+0.09 
-0.07 



0.89 
0.98 



1 25+°" 

1 9Q+0.07 
^■^^-0.06 
1 Q4+0.04 
^■•'^-0.03 



20.79!^™ 
21.37!2.45 

34.53!;:^« 



1 1 9+20 
13.9+|-^ 
16.8!j;^ 



13.47+1;^^ 
29.08!l:52 



17 55+3-0° 

17 41+1-61 
^ '-^^-0.99 
91 9Q+1.09 
^^■^^^-1.02 



9.71 



-1.58 
-1.29 

10 59+0-69 



+ 1.30 



10.97 



11.05!2-44 

10.78![:°^ 
17.07^^^ 



13 67+ '-°3 


ii.oo!2-^« 





124 


1 ^ Q4+0-64 
^•^■^ -0.52 


10.68!|;^1 





133 


13 29+°-2^ 


10.30!°j« 





129 



8.50!°:^« 
8.30!°i2 



8.21 



+0.77 



7 07+2.50 
'■-"-1.57 





115 


c 40 +0.46 
-"■^^-0.47 





151 


7 30+2-'° 

'■-^"-1.51 





113 



10 26+°-^'' 


6 62+'-''° 

"■""^-1.42 





155 


10 24+°-''^ 


7 49+1-29 
'■^^-1.00 





138 


11 05+°-^'' 


8 14+°-''^ 

"■^^-0.59 





136 



+0.020 

-0.017 
+0.009 
-0.012 
+0.005 
-0.005 



+0.022 
-0.023 
+0.008 
-0.008 
+0.018 
-0.018 



+0.026 
-0.023 
+0.015 
-0.018 
+0.007 
-0.007 
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u— V Radius (kX) 



Fig. 1.— SZ profiles for A1835 (left), CL1226 (middle), and A1914 (right), plotted as a function 
of M,v-radius ( Vm^ + v^). The panels show the real component of the visibilities over the full 
range of the fit models, using a broken axis to capture the model predictions as the M,v-radius 
approaches zero kA. Using Eq.[T8l the visibilities were rescaled to Y{u, v) d\, which removes both 
the frequency and redshift dependence of the cluster visibilities. The black points with error bars 
(1-cr) are the binned 30-GHz SZA visibility data. Models for the compact sources have been 
subtracted from the visibility data before rescaling. The blue solid line is a high likelihood NOV 
model fit, while the red dashed line is a similarly-chosen fit of the isothermal /3-model. For the 
available data points, both cluster models fit equally well. However, note that as the M,v-radius 
approaches zero kA (corresponding to large spatial scales on the sky) — where there are no data to 
constrain the models — the jS-model predicts much more flux than the NOV model. The middle three 
(green) data points for CL1226 are taken from 90-GHz SZA observations. 
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Radius (arcsec) 

Fig. 2. — X-ray surface brightness profiles for A1835 (left), CL1226 (middle), and A1914 (right). 
The vertical dashed line denotes the 100 kpc core cut. In each panel, the blue, solid line is the 
surface brightness computed using a high-likelihood fit of the N07-I-SVM profiles (analogous to 
the SZ models plotted in Fig.[T]), while the red, dot-dashed line is that from an isothermal yS-model 
fit. Both fit lines include the X-ray background that was simultaneously fit with the cluster emission 
model. The black squares are the annularly-binned X-ray data, where the widths of the bins are 
denoted by horizonal error bars. The vertical error bars are the l-cr errors on the measurements. 
Arrows indicate r25oo and rsoo derived from the N07-I-SVM profiles (see Tables |4]& [5]). 
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Fig. 3. — Three-dimensional ICM radial profiles derived from the joint analysis of SZA+Chandra 
X-ray surface brightness data for A1835, CL1226, and A1914 (from left to right). From top to bot- 
tom: the electron pressure Peir), the electron density n^ir), the derived electron temperature T^ir), 
profiles as a function of the cluster-centric radius. The lines show the median deprojected quantity 
derived from data using the N07-I-SVM (solid lines, vertically-hatched regions) and the isothermal 
jS-model (dot-dashed lines, horizontally-hatched regions). The derived electron temperature pro- 
files are compared to the spectroscopically-determined radial temperature profiles (black dashed 
lines, slanted hatching) obtained according to methods presented in M07/M08. The hatched re- 
gion indicate the 68% confidence on parameters derived from each model. The arrows denote the 
median values of r25oo and rsoo from the N07-I-SVM fits. 
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Fig. 4. — Ycy\ computed within 6' (~ 1 .4 Mpc, which is ^ r^oo for this cluster) for fits to SZA 
observation of A1835. The bold, black contours contain 68% and 95% of the accepted iterations 
on the jointly-fit Chandra + SZA data, while the thinner, blue contours are those for fits to SZA 
data alone. 
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X 10"^ 




Normalized Parameter 

Fig. 5. — 1-D histograms of Mtot (cyan/light gray region with dashed outline) and Fcyi (vertically- 
hatched region with solid outline), normalized by their respective median values, derived from 
the joint fits of the NOV+SVM profiles to A1835. Both Mjot and Ycyi are computed within a fixed 
radius of = 6'. The derived Ycyi, which scales directly with integrated SZ flux, has a more tightly 
constrained and centrally peaked distribution than that of Mtot, as Mtot is sensitive to the change in 
slope in the pressure profile (dP/dr). 




Fig. 6. — Total mass profiles derived from the joint analysis of SZA+Chandra X-ray surface 
brightness data for A1835, CL1226, and A1914 (from left to right). The line types are the same as 
those in Fig. |3l 



